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ABSTRACT 

The CO-H2 conversion factor {Xco \ otherwise known as the X-factor) is observed to be re- 
markably constant in the Milky Way and in the Local Group (aside from the SMC). To date, 
our understanding of why Xqq should be so constant remains poor. Using a combination of 
extremely high resolution 1 pc) galaxy evolution simulations and molecular line radiative 
transfer calculations, we suggest that Xqq displays a narrow range of values in the Galaxy 
due to the fact that molecular clouds share very similar physical properties. In our models, this 
is itself a consequence of stellar feedback competing against gravitational collapse. GMCs 
whose lifetimes are regulated by radiative feedback show a narrow range of surface densities, 
temperatures and velocity dispersions with values comparable to those seen in the Milky Way. 
As a result, the X-factors from these clouds show reasonable correspondence with observed 
data from the Local Group, and a relatively narrow range. On the other hand, feedback-free 
clouds collapse to surface densities that are larger than those seen in the Galaxy, and hence re- 
sult in X-factors that are systematically too large compared to the Milky Way's. We conclude 
that radiative feedback within GMCs can generate cloud properties similar to those observed 
in the Galaxy, and hence a roughly constant Milky Way X-factor in normal, quiescent clouds. 
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1 INTRODUCTION 

Historically, determining H2 gas masses in the interstellar medium 
(ISM) of galaxies has relied on the usage of tracer molecules. This 
owes to the fact that H2 requires temperatures ~ 500 K to excite 
the rotational lines, and is thus a poor tracer of cold (~ 10 — lOOK) 
giant molecular clouds (GMCs). Carbon Monoxide (^^CO; here- 
after, CO) is the second most abundant molecule, and has strong 
lines in a readily accessible window of atmospheric transmission. 
For this reason, CO is the most commonly employed tracer of H2. 
However, utilising CO as a measure of H2 does not come without 
uncertainty. 

At the heart of converting CO line fluxes to H2 masses is 
the CO-H2 conversion factor. The conversion factor relates ei- 
ther H2 column densities to velocity-integrated CO line intensity 
(-^co ; alternatively, the X-factor), or H2 gas mass to CO line lu- 
minosity (qco)- The two are related via Xco (cm~'^/K-km s~^) = 
6.3 X 10^^ X QCO (Mq pc"^ (K-km s"^)"^). Uncertainties in CO 
abundances, H2 gas fractions, and radiative transfer all complicate 
our understanding of Xco • 

In principle, Xco can be empirically calibrated with an in- 
dependent measure of H2 gas masses. Efforts along these lines 
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have used CO line widths (combined with an assumption regard- 
ing the dynamical state of the GMC), dust mass measurements 
(combined with an assumed dust-to-gas ratio), or gamma-ray ob- 
servations from GMCs to determine the CO-H2 conversion factor 
(Larson" 1981'; 'Solomon et al." 1987; 'Dickman ' 1975'; 'de Vries et alj 
1987; Pineda etal. 2008; Leroy et al. 2011; Magdis et al. 20ll|; 
Bloemenetal. 1986; Strong & Mattoxl 1 19961 ; lAbdo et al.l i201ol ; 
iDelahave et aL.201 1) . These measurements have all arrived at the 
conclusion that the X-factor in Milky Way GMCs is remarkably 
constant, displaying a relatively narrow range of Xco ~ 2 — 4 x 
lO^^cm^^/K-km s"^ (qco ~ 3 - 6 Mq pc~^ (K-km s"^)"^). 
Beyond this, observations of GMCs in relatively normal galaxies 
within the Local Group (that is, excluding the Small Magallenic 
Cloud; lEerov et al]|201 ih have evidenced similar X-factors as in 
the Galaxy. 

This said, not all galaxies exhibit X-factors comparable to 
the relatively constant Local Group disc galaxy value. In particu- 
lar, heavily star-forming systems at low and high- 2; appear to have 
X-fa ctors roughly a factor 2-10 lower than the Milky Way average 
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Downes & Solomon|[T993 , iTacconi et al. I2OO8I ; iMeier et al.1 
Naravananll20I ih . whereas low-metallicity galaxies can have 



fe.a. IWilsoiJll995l;lArimoto et al.||l99e 


; Israel 


1997; 


Boselli et al. 


2OO2I; iBolatto et al.f 12008; .Lerov et al.l 


201 ll; 


Genzel et al.l |2012|; 


Schruba et al. 20121). Critical questions include: d) what is the ori- 



gin of the nearly-constant X-factor in nearby disks, (2) why is Xco 
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depressed in high gas-surface density environments, and (3) why is 
the conversion factor elevated in low metallicity galaxies? 

In recent years, there has been a flurry of interest from the- 
orists in understanding Xco on scales ranging from GMCs to 
cosmological simulations of galaxy formation. Models have made 
great headway in understanding the latter two questions. MHD cal- 
culations of evolving GMCs (Shetty et al. 201 1a, b), cosmological 
galax y formation calculations llFeldmann et al.ll2012l : ILagos et al.l 
l2012i) and hydrodyn amic galaxy evolution calcula tions coupled 
with radiative transfer lNaravanan et alj ( [2011bll2012h have all con- 
verged on a picture in which low-metallicity galaxies have large 
fractions of CO-dark molecular gas (due to photodissociation of 
CO in regions of low dust extinctior[^ , and hence large X po- 
Similarly, galaxy merge r models by ,Naravanan et alj ( 1201 id) and 
[Narayanan et alj ( 1201 2h have shown that starburst environments 
can force large gas temperatures and velocity dispersions which 
increase the CO intensity at a given H2 gas mass, thus reducing 
XcQ. 

This said, thus far no model that considers a full ensemble 
of clouds on galaxy-wide scales has been able to explain why the 
X-factor has such a narrow range of values in the Milky Way 
and nearby galaxies. To understand this requires knowledge of the 
physical state of GMCs on highly-resolved pc) scales, but sam- 
pled over the scales of entire galaxies. That is, one ideally should 
be able to super-resolve GMCs while capturing the effect of the 
larger galactic environment on cloud evolutiorQ Moreover, with- 
out explicit models for feedback, GMCs will inevitably collapse 
without limit (to arbitrarily high densities) and turn most of their 
mass into stars. This is in stark disagreement with observations in- 
dicating in efficient star formation and relatively short GMC life- 
times (e.g. lEva ns 1999: Krumholz & Taj 20071 : Ishirlev et alj2007l : 
iKennicutt & Evans.2012). 

Recently, iHopkins et aD (I20T ij) have implemented various 
forms of stellar feedback into idealised galaxy evolution simula- 
tions that have allowed for ~ pc-scale resolution on galaxy-wide 
scales. These simulations move beyond previous models that em- 
ploy subre solution assumptions governin g molecular cloud evo- 
lution (e.g. [Narayanan erSll201 Ibl 1201 2h and allow us to super- 
resolve GMCs on galaxy-wide scales. These simulations have been 
utilised by Hopkins et al. (2012b) to show that a model in which 
radiative feedback from massive stars dominates the life-cycle of 
molecular clouds successfully reproduces many observed physical 
properties and scaling relations (e.g. "Larson's Laws") of GMCs, 
including their surface densities, velocity dispersions, size distri- 
butions and mass spectra. Here, we employ these simulations to 
model a Milky Way-like disc galaxy and ask whether GMCs whose 
physical properties are governed by radiative feedback can explain 
X-factor properties (i.e. their constancy) of observed Local Group 
GMCs. 



^ This is the presu med origin for the large observed X-factors in the SMC 
iLerov et al.l20ril). 

^ In lNaravanan et al.l feOllbl) , we implemented a subresolution model for 
the velocity dispersions and surface densities of GMCs in idealised galaxy 
evolution models. While these models suggested that the origin of a roughly 
constant Galactic X-factor owed to a relatively narrow range of physical 
properties in the model GMCs, this conclusion was not entirely independent 
of subresolution assumptions. 



2 METHODS 

We simulate the hydrodynamic evolution of a Milky Way- 
like galaxy with a substantially modified version of GADGET-3, 
a smoothed-particle hydrodynamics (SPH) code (^Spring el et al.l 
2005: Springel 2005). The mai n code modificatio ns and model de- 
tails are described in Hopkins et al I ( l201lLi2012bl) and we describe 

only the important aspects here. 

W e initialise an exponential disc according to the IMo et al.l 
( Il998 h model within a live dark matter halo of mass 1 .6 x IO^^Mq 
with a Hernquist ( 199i)) density profile and concentration parame- 
ter c = 12. The baryon, bulge, disc and gas masses are initialised 
at (Mbar, mi,, rrid, nig) = (7.1, 1.5, 4.7, 0.9)x IO^^Mq, and scale- 
lengths (/id, ftj, zo)=(3.0, 6.0, 0.3) kpc. The simulations are run 
with ~ 3 X 10^ hydrodynamic particles, with a smoothing length 
of ~ Ipc. 

For the purposes of the SPH simulations, the gas is allowed 
to cool to ^ 100 K (though see below regarding further refine- 
ments on this temperature structure in post-processing), and the 
HI-H2 balance i s determined following the semi-analytic model of 
iKrumholz et al. 1(12008, ,,,2009). Stars form exclusively in H2 gas that 
is self-gravitating on the smallest resolved scales. Stars form at an 
instantaneous rate of p« = p/tg{p); because feedback can prevent 
further star formation once stars form, the average efficiency of star 
formation in dense gas is typically ~ 1%. 

The most important aspect of the models is a range of mecha- 
nisms by w hich stellar feed back can impact the ISM. Throughout, 
we assume a'Kroupa ('2002') stellar IMF, and utilise STARBURST99 
(Leitherer et al. 1999) for all stellar luminosity, mass return and su- 
pernova rate calculations as a function of stellar age and metallicity. 
The mechanisms of stellar feedback include: 

(i) Local Momentum Deposition: For the purposes of this 
work, the most important source of feedback is local momentum 
deposition by stellar radiation, mass return from stellar winds and 
supernova. At each timestep, the nearest density peak to a given 
gas particle is determined to represent a clump inside of a GMC. 
The total stellar radiation from all star particles inside a sphere de- 
fined by the distance between the gas particle and density peak are 
summed; the radiation from these stars is then used to determine 
the momentum flux. 

The momentum flux from radiation is given by Prad ~ (1 + 
Tm) X (L/c), where tir = EgasKra,. Sgas is the gas surface density 
calculated directly from the simulation, and Km. is approximated by 
KIR = 5{Z/Ze)g'^ cm^. 

We also include direct momentum injection from SNe and stellar 
winds, whose momentum deposition rates are directly calculated 
from STARBURST99 and injected to the gas within a smoothing 
length of the star particle. This source of feedback is typically sub- 
dominant compared to radiation pressure in galaxies as massive as 
a MW analog. 

(ii) Thermal Energy Input from Supernova and Stellar 
Wind Shock-Heating: For the supernovae (SNe), we tabulate Type 
I and II SNe rates from Mannucci et al. (2006) and STARBURST99, 
respectively, for all star particles and determine from a stochastic 
procedure if a SN occurred at each timestep. When a SN occurs, 
thermal energy is injected into the gas within a smoothing length 
of the star particle. For stellar winds, we inject the tabulated me- 
chanical power as a function of stellar age and metallicity to the 
gas within a smoothing length of the star. 

(iii) Photo-Heating of HH Regions The production rates of 
ionising photons from star particles is calculated and used to de- 
termine the extent of HII regions surrounding stars (allowing ap- 
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propriately for overlapping regions). The temperatures of these HII 
regions are heated to 10'* K if the gas falls below that threshold. 

(iv) Long-Range Radiation Pressure Photons that escape the 
local GMC (not accounted for in (i)), after being appropriately 
attenuated/absorbed, are propagated to large distances along di- 
rect rays, where local absorption is calculated by integrating over 
a frequency-dependent opacity (which also scales linearly with 
metallicity). The appropriate radiation pressure forces are then im- 
parted. 

Predicting the molecular emission requires calculating the cor- 
rect temperature (including full radiative transfer effects) of molec- 
ular gas below T — 100 K (the approximate floor imposed by the 
cooling tables in the simulations). This is prohibitively expensive 
to perform on-the-fly, so in this paper we do so in post-processing. 
But we do not expect changes in the thermal pressure at such low 
temperatures to have any dynamical effect on the simulations. We 
calculate the temperature utilising the methodology described in 
[Goldsmith (2001) and Rrumholz et al. (201 lb). The dominant heat- 
ing processes of the H2 gas are cosmic ray heating, the grain photo- 
electric effect, and dust-gas thermal exchange. The dominant cool- 
ing terms are cooling by CO, CII and dust. If the gas and dust are 
in thermal balance, then we have the following equations (where F 
represents heating terms, and A represents cooling terms): 



+ TcR — Alin, 
Tdust A(jus 



(1) 

(2) 



The equation is solved by simultaneously iterating on the tempera- 
tures of the gas and dust 

We assume a Galactic cosmic ray heating rate and a grain pho- 
toelectric eff ect proportional t o the lo cal FUV intensity. We refer 
the reader to iKrumholz et al.l ( 1201 lal) for the specific values em- 
ployed in the model. In short, however, the temperature can be 
thought of as density-dependent. At high densities (n > lO^cm^"^), 
dust and gas exchange energy efficiently, and the gas temperature 
rises to the dust temperature. At low densities (n < lO'^cm"'^), cos- 
mic rays dominate the heating. For a Milky Way cosmic ray flux, 
this corresponds to GMCs with temperatures ~ 8 — 10 K. Interme- 
diate densities have temperatures in between the dust temperature 
and cosmic ray-determined temperature. 

The dust temperature is calculated utilising SUNRISE, a pub- 
licly available Monte Carlo dust radiative transfer code (see 
Ijonsson et al.ll201ol: Jonsson & Primack''20I0. for code descrip- 
tions, as well as iHavward et alj 1201 k i2012ai ib| ) for further de- 
tails). W e utilise the simulation set-up described in lNaravanan et al] 
( l201Ibl) , and refer the reader there for more details. In practice, the 
bulk of the gas remains below density n = lO^cm^'^, the density 
at which dust-gas energy exchange becomes efficient. 

The line cooling te rm is calculated in each cell via a ID 
escape probability code dKrumholz & ThompsonI |2007|) . We as- 
sume a fractional carbon abundance of 1.5 x 10~*Z' , where Z' 
is the metallicity with respect to solar. The fraction of hydro- 
gen where the carbon is in the form of CO is well approximated 
from both semi- analytic dmifire et alJl2oTol) and numerical mod- 
els (|GioveL&M^^^20il|): 
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where Gg is the FUV intensity relative to the Solar neighbourhood. 
Physically, Equation[3]describes the photodissociation of CO from 
UV photons, and its ability to survive behind sufficient columns of 



dust. When fco > 0.5, we assume that CO dominates the line 
cooling; else, CII. 

With the physical and chemical state of the molecular gas 
known, we utilize TURTLEBEACH, a 3D non-local thermody- 
namic equilibrium adaptive mesh Monte Carlo line radiative trans- 
fer code to calculate t he velocity in tegrated CO line intensity 
(Narayanan etal. 200|, I2OO8I l20Ilal lbl). We refer the reader to 
iNaravanan et alj ( 1201 IW for the formal equations, and only sum- 
marise the relevant points here. 

CO line emission is set by the level populations. The source 
function for a given transition it ^ Z is given by: 



(4) 



(niBiu — UuBui) 

where A^^i , Bi^ and B^i are the Einstein rate coefficients, and n 
are the level populations. 

We first calculate the level populat ions within a given cell 
utilising the escape probability formalism (' Krumholz & ThompsonI 
2007). We emit model photons from each cell isotropically with 
emission frequency drawn from a Gaussian profile function. When 
the photon passes through a cell, it sees an opacity of: 



Q"'(gas) = ^^^(l>{iy)(niBin 
47r 



nuBui 



(5) 



where u is the transition frequency, and <^(i^) is the line profile 
function that takes into account the effects of line of sight velocity 
offsets in the opacity. 

Once the model photons have all been emitte d, the level popu- 
lations are updated by assuming detailed balance dNaravanan et al.l 
l20Ilbh . The coUisional rate coefficients are taken from the Leiden 
Atomic and Molecular Database ( Schoieretalj2005h . This process 
is iterated upon until the level populations are converged to < 1% 
across all cells. 

GMCs within the model are identified via a friends of friends 
finder with a linking length of 20% of the mean cell size. Tests have 
shown that the results are not substantially sensitive to this choice. 



3 WHY IS THE X-FACTOR CONSTANT? 

To first order, Xco can be thought of as the column density of 
GMCs divided by the product of their temperature and velocity dis- 
persion. Formally, Xco= A'h2/Wco, where Wco is the velocity- 
integrated CO intensity. When the gas is in local thermodynamic 
equilibrium (as CO J=I-0 almost always is), the amplitude of the 
emission line is proportional to the gas kinetic temperature Tr. 
Similarly, because CO (J=l-0) is typically optically thick within 
GMCs, increasing the velocity dispersion of the gas increases the 
emergent CO intensity. As a result, Wco increases for both increas- 
ing kinetic temperature, as well as increasing velocity dispersion. 
So, to ask why the X-factor in MW GMCs is nearly constant is to 
ask why observed gas temperatures, velocity dispersions, and sur- 
face densities have a narrow distribution of values. 



3.1 Tlie Physical Properties of GMCs in Galaxy Discs 

In Figure [T] we show the mass-weighted distributions of GMC 
temperatures, velocity dispersions, surface densities, and A-factors 
for our fiducial MW model at a randomly chosen time snapshot 
(though the results are consistent for the bulk of the galaxy's evo- 
lution). As we will see, the physical properties of the GMCs are 
generally determined by the radiative feedback that eventually dis- 
rupts the GMC. 
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Figure 1. Mass-weighted distributions over all GMCs of gas surface density, velocity dispersion, temperature, and Jf-factor for our fiducial model (solid 
black line), and a feedback-free model (dotted line). Feedback-free GMCs collapse to higher surface densities than models that include feedback, driving more 
power into large Xqq values. 



Quiescent GMCs in the MW can be thought of as roughly 
isothermal, with temperatures ~ 8 — 10 K. The average densities 
of the GMCs are relatively low, with n ^ 10 - 100 cm~^. At 
these densities, the gas is energetically decoupled from the dust, 
and cosmic rays act as the primary heating source. For a Galactic 
cosmic ray flux (which we assume), this equates to a minimum gas 
temperature of roughly ~ 8 — 10 K. 

The GMC surface densities display a range of values cen- 
tred around ^ 100 Mq pc~^. We calculate the surface density as 
Eh2 = M/(tv X i?c), where the radius of the GMC (Rc) is half 
the average of the maximum length in three orthogonal directions. 
While this is method is not free of geometric effects, it is a rea- 
sonable approximation to the methods used in observations. These 
surface densities are comparable to collapse conditions for gas that 
forms GMCs. In a model where the GMC lifetime is regulated by 
radiative feedback, GMCs do not collapse indefinitely. Once GMCs 
reach surface densities near 100 Mq pc^^, feedback from star for- 
mation disperses the cloud. These surface densities are not far from 
the average value of the disc, as the model GMC only lives a few 
X 10^ yr (e.g. a few free fall times). Beyond this, these surface den- 
sities are comparable to those seen in Galactic GMCs, which dis- 
play a relatively n arrow rang e (Larson 1981; Solornon et al.i, 1987. ; 
ikever et alj2009l though see lLombardi et alj ( l2010h ). 

It is worth noting here that the X-factor is implicitly depen- 



dent on the volumetric density of the GMCs residing within a range 
of n ~ 50 — 10* cm~^. At larger densities, the gas and dust 
exchange energy efficie ntly, and the gas t emperature rises to tha t 
of the dust temperature dColdsm ith 2001; Krumholz et al.|[20I Ibl) . 
driving it to larger values than the roughly ~ 10 K value seen when 
cosmic rays dominate the heating. At lower densities, the CO may 
not be in LTE (depending on the degree of line radiative trapping). 
In this regime, the peak CO intensity will no longer scale with the 
kinetic temperature of the gas. As in the case with the GMC surface 
densities, radiative feedback sup presses the f ormation of excessive 
amounts of dense gas (.Hopkins et al.ll2012al) . The median density 
of GMCs is ~ 100 cm~^, and the distribution within the galaxy is 
roughly lognormal in shape. 

The GMCs in these simula tions are consistent w ith being 
marginally gravitationally bound dHopkins et al.ll20I2bl) , and have 
velocity dispersions ranging from a few to ~ 20 km s~^. Disrup- 
tion of the GMC by radiative feedback keeps the GMC from be- 
coming too strongly self-gravitating, and hence limits the velocity 
dispersions. Comparing t hese to the typical ve locity dispersion of 
GMCs seen in the Galaxy jSolomon et al.l 19871) . the range of mod- 
eled GMC velocity dispersions is reasonable. 



Given a Milky Way cosmic ray flux. 
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Figure 2. X-factors of model GMCs in our fiducial model galaxy as a 
function of galactocentric distance. Black circles represent the individual 
GMCs, while the red filled circles denote the median of all GMCs within 
bins of ~ 3 kpc (with the dispersion within the bin shown as error bars). 
The blue shaded region shows the rough range of values for the Milky Way. 
GMCs toward the centre of the galaxy tend to have lower X-factors due 
to systematically larger velocity dispersions, though the trend is relatively 
weak, and there is significant dispersion. 
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Figure 3. Virial mass versus CO luminosity for observed GMCs (red 
stars), and modeled clouds (black filled circles). Left plot shows our fidu- 
cial galaxy, and right plot a feedback-free model. While the clouds in our 
fiducial galaxy show excellent correspondence with observed GMCs, the 
feedback-free clouds systematically has X-factors that ai'e too lai'ge, driv- 
ing the simulated Afvir — ^CO relation a factor of a few higher than the 
observed points. 



3.2 X-factor properties in Galaxy discs 

We can now see why the X-factor is relatively constant in Milky 
Way-like disc galaxies. If the evolution of GMCs is largely gov- 
erned by radiative feedback, their surface densities and velocity 
dispersions display a relatively narrow range of values comparable 
to measurements of Galactic GMCs. If, beyond this, the tempera- 
tures of GMCs are nearly isothermal at ~ 10 K, as expected for 
clouds where the dominant heating source is cosmic rays with a 
flux comparable to the Galaxy's, then the CO-H2 conversion factor 
will display a relatively narrow range of values centred around 2-4 
X 10^" cm^^/K-km s^^. We see this in the fourth panel of Fig- 
ure [T] 

In Figure [2] we show the relationship between Xco in our 
model GMCs and their distance from the centre of the galaxy. The 
black points denote the individual GMCs, while the red circles 
and error bars denote the median and dispersion within roughly 
3 kpc bins. The blue shaded region shows the average Galactic 
range for the X-factor. Generally, GMCs toward the centre of the 
model Milky Way analog systematically have lower X-factors than 
clouds at larger distances from the galactic centre. Influenced by a 
large stellar potential, as well as other gas, GMCs toward the cen- 
tre of the galaxy tend to have larger velocity dispersions than field 
GMCs, and tend to be unvirialised. There is some tentative observa- 
tional evidence that Xco values in GMCs decrease from t he Ga lac- 
tic mean toward the centre of the Milky Way foka et alll 19981) . as 
well as in other nearby galaxies (Sandstrom et al. in prep.). This 
said, there is significant dispersion in observed trends of Xco with 
galactic radius, and some observations show nearly no depression 
at all (less than a factor of 2) toward galactic nuclei (Donovan 
Meyer et al. in prep.). 

The X-factor from our model GMCs compares well to the 
observed population of GMCs within the Galaxy. In the left panel 
of Figure [3] we show the virial masses of observed GMCs in 



the Galaxy and NGC 69 46 against their CO lumino sities from 
ISolomon et alj l ll987h and lDonovan Meyer et 'Zl( l2012h . We addi- 
tionally show the Mvir — Leo relation for our model galaxies. We 
defer discussion of the right panel for § 13.31 The normalisation of 
the A/vir — Leo relation in the model GMCs (black solid points), 
which betrays the conversion of CO to H2 gas mass, corresponds 
reasonably well with the observed data (red stars). 



3.3 Feedback-Free Models 

In order to highlight the role of feedback in setting the physical 
properties of the model GMCs, it is worth considering the proper- 
ties of a feedback-free model. At this point, we now highlight the 
dotted red line in Figure[T] which denotes a model run with exactly 
the same initial conditions as our fiducial model, though with no 
forms of feedback included. When GMCs first form, their proper- 
ties are not especially different in models with and without feed- 
back, as expected fro m general mod els of gravitational instabil- 
ity and fragmentation I^Ho Dkinsll20Ill) . But GMCs in models with- 
out feedback proceed to collapse one-dime nsionally, developing a 
pancake-shaped geometry and spinning up ( Ho pkins et al.ll2012bl) . 
While the general surface density distribution is not terribly dissim- 
ilar from our fiducial model, feedback-free GMCs develop a large 
tail of very high-density gas, and show more power toward large 
Eh2. The velocity dispersion distribution additionally increases 
moderately (due to a spinning up of the contracting GMC), as well 
as the temperature (due to more gas at or above the gas-dust cou- 
pling density of ~ 10*cm~^), but these increases are more modest 
than the increase in the GMC surface densities. 

The net result of feedback-free GMCs is more power toward 
large X-factors, and typical X-factors a factor ~ 2 — 3 larger 
than typical MW GMCs. As an example, in the right panel of Fig- 
ure [3] we plot the Afvir — Leo relation for observed GMCs and 
model GMCs in a feedback-free model. Compared to our standard 
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model, which shows excellent correspondence with observed data, 
the feedback-free model lies a fa ctor of a few above the M yir — Leo 
relation. Beyond this, as noted bv H opkins et alj ( l2012bh . feedback- 
free models exhibit a host of problems, including star formation 
rates well above the observed Kennicutt-Schmidt relation, gas den- 
sity distributions (and HCN/CO line ratios) above those observed 
jHopkins et alj|20^12al) . and abnormally low virial parameters. 

3.4 Summary 

Utilising a combination of high-resolution galaxy evolution simu- 
lations that super-resolve giant molecular clouds and 3D molecular 
line radiative transfer calculations, we investigate why the CO-H2 
conversion factor {Xco) is observed to be roughly constant in the 
Milky Way and Local Group galaxies (aside from the SMC). 

Our main result is that Xco is found to be nearly constant be- 
cause all GMCs in our model Milky Way have similar physical con- 
ditions to one another. In particular, Xco is determined principally 
by GMC surface densities, temperatures, and velocity dispersions. 

• The model GMCs all have similar surface densities at val- 
ues near ~ 100 Mq pc^^. Models that include radiative feedback 
limit GMCs from achieving surface densities much larger than this 
value due to cloud dispersal. 

• The temperatures of GMCs are dominated by cosmic ray heat- 
ing. Given a Milky Way cosmic ray flux, this results in nearly 
isothermal GMCs with temperatures ~ 8 — 10 K. 

• The GMCs in our model are consistent with being marginally 
bound, resulting in a narrow velocity dispersion range of 5 — 20 
km s"'^. 

These GMCs display a relatively narrow range of physical 
properties, and compare well with those observed in the Milky Way. 
As a result, the CO-H2 conversion factor in these clouds addition- 
ally shows a narrow range of values. Feedback is a necessary ele- 
ment in the model to control the cloud surface densities, and hence 
limit the observed range of Xco- Feedback-free models collapse 
to large surface densities, and hence show excessive power to large 
Xco values. 
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